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The spin-1/2 quantum anisotropic XY spin chain in a transverse random magnetic field parallel 
to the z axis is numerically studied by means of the density-matrix renormalization group. The 
dependence of the spontaneous magnetization and the energy gap on both the strength of the random 
magnetic field and the factor of anisotropy is determined. The critical line for the order-disorder 
phase transition is obtained and the resulting phase diagram is drawn. Our results are compatible 
with the fact that models with different factors of anisotropy fall within the same universality class 
as the quantum Ising model in a transverse random field. 



PACS numbers: 75.10.Jm, 75.40. Mg, 05.50.+q 
I. INTRODUCTION 

The effect of disorder on quantum systems is more in- 
triguing than its classical counterpart, especially when 
the random noise does not directly affect the interactions 
in the original system, but is coupled to other degrees 
of freedom, which have nontrivial commutation relations 
with the degrees of freedom appearing in the original 
Hamiltonian. A simple example within this class of sys- 
tems is provided by the transverse randomness in an easy- 
axis or easy-plane quantum spin chain. Indeed, when 
the exchange interactions involve the x and y spin com- 
ponents, a random magnetic field coupled to the z spin 
component produces dramatic effects due to the nontriv- 
iality of the angular-momentum algebra. 

In this paper we analyze the one-dimensional spin-1/2 
quantum anisotropic XY model in a transverse random 
magnetic field, which is defined by the Hamiltonian 
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are the spin-1/2 operators acting on the Hilbert space 
spanned by the vectors . . . , tpe, . . . , i/jn), {4>l} are 
two-component spinors, <r° is the 2x2 identity matrix, 
{o- a } are the 2x2 Pauli matrices, ® is the symbol for 
tensor product, J is the coupling constant for the x spin 
components, 7 is the factor of anisotropy, 7 J is the cou- 
pling constant for the y spin components, and {hi} are 
the on-site values of a random magnetic field directed 
along the z axis, which are uniformly distributed in the 



interval —ho<hi<ho. In the following we call the 
strength of the random magnetic field. 

Due to the conventions adopted in Eq. (|l|) a positive 
(negative) coupling constant stands for a ferromagnetic 
(antiferromagnetic) coupling. Correspondingly the exter- 
nal field h x directed along the x axis, which we include 
in our model to force spontaneous symmetry breaking 
along the easy axis of magnetization, is uniform (stag- 
gered) when J is ferromagnetic (antiferromagnetic). 

In principle the coupling constant J and the factor 
of anisotropy 7 have arbitrary signs, and the only im- 
plicit limitation in Eq. (Jl]) is — 1 < 7 < 1, which 
amounts to defining the x and y axes according to 
the relative strength of the exchange interaction. How- 
ever the local gauge transformations Sf — > (sgnJ)^S 1 ^, 
Sf — > (sgnJsgn7)^S']', and Sf — ► (sgaj) Sf make all the 
coupling constants ferromagnetic and the external field 
h x uniform, while preserving the commutation relations 
among the spin operators, i.e. the orientation of the local 
reference frames is preserved by the transformations. We 
take then the units of energy such that J = 1 , and only 
consider the case < 7 < 1. We also assume that h x > 0. 
It is important to observe that the gauge transforma- 
tions change the random field as he — » (sgnj) h^ which 
is again uniformly distributed in the interval [— ho, ho]. 
However, unlike the Ising limit (7 = 0), at finite 7 there is 
no freedom left to make all the {hi} positive (cf. Ref. [l]), 
because any further local gauge transformation will not 
keep fixed both the sign of the coupling constants and 
the orientation of the local reference frames. 

The outline of the paper is the following. In Sec. II 
we briefly recapitulate some previously established re- 
sults in the two limiting cases: the XY model in the 
absence of randomness and the quantum Ising model in 
a transverse random field. The main scope of Sec. II 
is to provide a framework to which to refer in the dis- 
cussion and the interpretation of the results obtained in 
this paper. In Sec. Ill we discuss some technical aspects 
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which are fundamental to obtain well-defined numerical 
results, and their relation to the physical properties of 
the model. This discussion is essential to allow for the 
reproducibility of our results, and to shed light on the 
physics underlying the DMRG algorithm. In Sec. IV we 
determine the critical properties and the phase diagram 
of the model (jl|) in the ho vs 7 plane. A summary of the 
principal results and some concluding remarks are found 
in Sec. V. 



II. SOME LIMITING CASES 

The model (§) with a generic 7 is integrable in the ab- 
sence ofjijandomness and external magnetic field = 0, 
h x = Op by fermionization of the spin operators .□ The 
particle-hole excitation spectrum of the corresponding 
fermionic system shows, in particular, that the ground 
state is separated from the first excited state by an en- 
ergy garjH 
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(in units of J) . In the thermodynamic limit the spectrum 
above the gap is continuous and the ground-state energy 
per site is 
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A pictorial representation of the physics of the XY model 
is provided by the spin wave (SW) theory, which is qual- 
itatively correct (in a sense that will be clearer in the fol- 
lowing), except for 7 very close to 1. Expressing the spin 
operators in terms of auxiliary boson operators at,(x[ 
according to Sf = \ — a\at, S\ = \(al + at), and re- 
taining only bilinear terms, the Hamiltonian (Q) is cast in 
a diagonal form by a Fourier transform from real-space 
to fc-space representation, and by a subsequent canonical 
Bogoliubov transformation. The resulting SW excitation 
spectrum is Ek = y/l — 7 cos k, and the ground-state en- 
ergy per site in the thermodynamic limit is 
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I.e., the gap in the excitation spectrum is G^w = yT — 7, 
which coincides with (Q) only in the Ising limit 7 = 
and in the isotropic limit 7 = 1, and is not even pertur- 
batively correct at small 7. Nonetheless the discrepan- 
cies are mostly quantitative and the SW theory may be 
adopted as a reference theory to discuss the physics of 
the XY model. In particular, at 7 > 0, the SW Hamilto- 
nian includes a contribution from the zero-point motion 
of elementary excitations. The reduction of the sponta- 
neous magnetization with respect to the saturation value 
Mo.max = \ in the Ising limit (7 = 0) is understood 



within SW theory as the effect of such quantum fluctua- 
tions for 7 > 0. Indeed the SW result for the spontaneous 
magnetization is 
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This result is meaningless at 7 = 1 due to the divergence 
of the integral in the r.h.s., a fact which is used as an 
evidence against the existence of a spontaneous magne- 
tization in the isotropic XY model. However, due to the 
logarithmic nature of the divergence, it turns out that 
the reduction of the magnetization associated with SW 
zero-point motion is an underestimate of the reduction 
in the exact ground state (see Eq. (^|) and text below it), 
except for 7 very close to 1. Thus in the following, for 
the sake of definiteness, we adopt the SW language and 
we conventionally refer to the effects of the anisotropy as 
resulting from quantum fluctuations with respect to an 
Ising-like ground state, as it would be in a more refined 
(i.e. beyond SW) field theoretical approach to the XY 
model. 

As a consequence of the gap in the excitation spectrum 
the quantum fluctuations are massive. The spins are or- 
dered and the system is ferromagnetic (FM) for 7 < 1. In 
the isotropic limit 7=1 (the so-called XX model), the 
gap closes and the spontaneous magnetization is driven 
to zero, by massless phase fluctuations. The system is 
thus paramagnetic (PM). It is worth mentioning that in 
the isotropic limit the model is equivalent to a model of 
hard-core bosons, and the absence of spontaneous magne- 
tization in the magnetic system is related to the absence 
of Bose-Einstein condensation in the corresponding one- 
dimensional interacting boson system. 

The issue addressed in this paper is the evolution of the 
magnetic ground state for 7 < 1 in the presence of the 
transverse randomness introduced by a random magnetic 
field along the z axis. In particular we want to study the 
interference between the critical behavior driven by 7 — > 
1 and the critical behavior controlled by the random-field 
ho. Indeed, in the extreme anisotropic case 7 = Ocffl thg 
system undergoes a FM-PM quantum phase transition.Q 
The magnetization at the critical point ho = h c = e/2 = 
1.3529 ... is a nonanalytical function of the external field 
h$ M (h = h c , h x ) ~ [In (D h /h x )f-\ where cf> = (1 + 
VS)/2 = 1.618 ... is the golden mean and Dh is a scale 
factor. The spontaneous magnetization Mq = M(h x — > 
0) in the ordered phase decreases with increasing ho, and 
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for ho — » h~ . Contrary to ordinary continuous phase 
transitions, where the nonanal*jticity is limited to the 
critical point, a Griffiths' regionB exists around the criti- 
cal point, where M(h x ) is not analytical. In the weakly 
ordered Griffiths' region (ho < h c ), for small h x , 



M(h x )-Mo~K[\n(D'Jh x )X 



(6) 



2 



where Mq is the spontaneous magnetization, the ex- 
ponent k cannot be determined within renormalization 
group, D' h is a nonuniversal scale factor, and the expo- 
nent a depends continuously on the distance from criti- 
cality 5 ~ h c — As far as a < 1, the spin susceptibility 
diverges as h x — ► 0, i.e. the magnetization curves inter- 
cept the axis h x = with an infinite slope, even when 
Mq is finite. 

Previous results of density-matrix Eenormalization- 
group (DMRG) calculations at » Otl were in good 
agreement botL-adth the theoryoQ and other numeri- 
cal calculations £3 Below we extend these results to the 
entire region ho > 0, and < 7 < 1, and determine 
the full phase diagram and the critical properties of the 
model (|l|). We are interested in understanding the role 
of quantum fluctuations (in the sense clarified above) in 
the Griffiths' region. The issue is particularly intriguing 
for 7 — > 1 , when the critical behavior of the nonrandom 
system interferes with the critical behavior controlled by 
the random field. This interference is directly evident in 
those quantities which vanish at criticality, such as the 
spontaneous magnetization and the gap in the excitation 
spectrum. 

We point out that, due to the iterative nature of 
the infinite-size DMRG algorithm, the system size is in- 
creased while the irrelevant states are progressively trun- 
cated away. The discrepancy between the infinite-size 
limit of the algorithm and the thermodynamic limit of 
the physical system depends_.on the size of the error in- 
troduced by the truncation.EJ Thus the technical aspects 
related to the definition of a proper numerical procedure 
are deeply entangled with the physical properties of the 
system. In particular it will be clear that, in order to 
drive the system towards the phase transition, it is not 
only necessary to reach the infinite-size limit, in the pres- 
ence of an infinitesimal external field, but it is also impor- 
tant to control and determine the effect of the truncation 
of the Hilbert space. The full procedure explained below 
is essential to reproduce our results and obtain the cor- 
rect limiting behaviors. 

III. ALGORITHM AND TECHNICALITIES 

We perform our calculations using—the infinite-size 
DMRG algorithm introduced by WhiteO with some mi- 
nor modifications described in Ref. [l]. Accordingly, we 
adopted open boundary conditions, which produce more 
accurate results within DMRG Jill 

At a fixed number s of states kept in the truncation 
of the density matrix and after a (sufficiently large) fixed 
number N-rq of DMRG steps, when the chain consists 
of N = 2Ar,g + 2 sites, the magnetization at a given 
(small) uniform magnetic field h x decreases both with 
the strength of the random field ho and the factor of 
anisotropy 7, except for a region close to 7 = 1 where 
the magnetization appears to be a nonmonotonic func- 



tion of ho (see Fig. [T]). The strong dependence of the 
magnetization on the external field around the critical 
point is essentially the result of a divergent susceptibil- 
ity in the Griffiths' region. Even for a field as small as 
h x = 10~ 5 , the transition point is masked by the strongly 
nonanalytical behavior, and is only hinted by the change 
in the curvature of the M vs ho curves. 




FIG. 1. The dependence of the magnetization on the 
random magnetic field for different values of the factor of 
anisotropy 7. The system size is 302 sites, the uniform mag- 
netic field h x — 10~ 5 , and eight states are kept at each DMRG 
step. The large value of the magnetization in the disordered 
region is essentially due to the anomalous response to the 
external field h x in the Griffiths' region. 

Moreover, it is evident that particular care has to be 
exercised close to the isotropic limit 7 = 1 , to extract the 
correct behavior. Indeed, quantum fluctuations suppress 
the spontaneous magnetization at 7 = 1 even in the ab- 
sence of randomness. Thus we expect Mo to be zero for 
all values of ho- This indicates that a proper procedure 
has to be defined, which is reliable and unique for all 7. 
We carefully address this issue and the results are given 
below. 

To determine the spontaneous magnetization a (small) 
uniform magnetic field h x in the x direction is applied to 
lift the degeneracy in the ground state, which is doubly 
degenerate at h x = (this degeneracy occurs for large 
system sizes for 7 > and for any size at 7 = 0). Usually, 
for small values of 7, any h x is enough to truncate the 
undesired degenerate state after the initial DMRG steps. 
However, the closer one gets to 7 = 1, the larger this 
uniform field should be, and more than 150 DMRG steps 
could be needed to isolate the state with broken sym- 
metry. Further, increased accuracy of calculations (i.e. 
a larger number of states kept) is also needed close to 
7 = 1. We have checked that the difference between the 
exact ground-state energy £o hco , given by Eq. (||), and the 
calculated one £(!p MRG for a fixed number of states kept 
in the truncation, procedure, s = 8, closely follows the 
truncation erroiEl 1 - P s = 1 - Y%n=i @ m = E m » 8m, 
i.e. the sum of the eigenvalues of the density matrix 
which correspond to the truncated states. This error in- 
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creases with increasing 7, from 10 -13 at 7 = 0, to 10~ 4 
at 7 = 1. This phenomenon is deeply related to the in- 
creasing complexity of the ground state, as the role of 
quantum fluctuations is emphasized. As a consequence, 
the spontaneous magnetization has a dependence on the 
parameters s and -/Vrg which becomes more and more 
important as 7 — * 1. 

To deal with this dependence and obtain the correct 
spontaneous magnetization Mq = M(h x — > 0) for all 
7, the following steps are needed: (i) We calculate the 
magnetization along the easy axis 

(N KG ,h x , h ) = (2N RG + 2)- 1 x 

(h x ,h ;N RG ,s,j\S x \h x ,h ;N RG ,sJ), (7) 

where ( • ) is the quantum expectation value in the ground 
state and S x is the properly truncated total x-spin oper- 
ator, at each DMRG step N RG , for a given uniform mag- 
netic field h x and a random magnetic field of strength 
ho, while s states are kept during the DMRG proce- 
dure. The index j labels a particular realization of the 
random field. Then (ii) we average this value over N c 

realizations of the random field, M (Nji G ,h x ,ho) = 
iV" 1 J2f=i Mj(NBG)h,c, ho), where a suitable value for 
N c depends on both 7 and ho- In most cases an aver- 
age over N c = 500 configurations gives satisfactory re- 
sults, though sometimes more than 2000 configurations 
are needed (for 7 ~ 1 and close to criticality). 

Due to the effects of the edges of the chain, M always 
increases with increasing lattice size, but the increase is 
very slow above some typical chain size (N RG = 50), and 

may be fitted as M (-/Vrg , h x , ho ) = AT '{ 00 , h x , ho) — 
( s \ 

A/JVrg, where A and M (00, h x ,ho) are fitting param- 
eters, and so (iii) we extrapolate M to the infinite-size 

limit Nr G — > 00 and obtain M * (00, h x , ho). 

The larger the number of states s kept is, during renor- 
malization, the better the accuracy obtained is, so the 
magnetization slightly decreases with s (a qualitative ar- 
gument is: the smaller the number of states kept is, the 
more influential the state with all spins up is, the larger 
M is). Thus (iv) we make an extrapolation to s = 00, til 
using the formula 

M <S \oa,h x ,h ) = M i ™\oc,h x ,h )+Bx (^J , (8) 

which fits the data quite well (except for small oscilla- 
tions at larger 1/s), and we obtain M °° (00, h x ,ho), (see 
Fig. ||). A more accurate formula is slightly different. On 
physical grounds, indeed, one expects that the magneti- 

zation M(l/s) = M (oo,h x ,ho) consists of two parts: 
M(l/s) = M reg (l/s) + Mnoise, where the regular part 
M reg (l/s) has a Taylor expansion close to 1/s = and 
the noise part M no i se is due to nonuniform changes in 
the ground-state properties and only exists at large 1/s. 



As the field h x is increased, the dependence on 1/s be- 
comes much weaker (in strong fields the ground state is 
fully magnetized and so s = 1 gives the exact solution), 
which means that for large fields one expects M(l/s) to 
be a noiseless horizontal line, which can be imagined as 
resulting from a successive vanishing of the terms in the 
Taylor expansion of M rog (l/s). Indeed, a Taylor expan- 
sion up to fifth order shows that only the fourth order 
1/s term has any influence for the h x = 4 x 10~ 3 data 
(Fig. ||) , though both the linear and square 1/s terms are 
present for h x = 1 x 10 -4 . The fact that no more than 
two nonconstant terms are ever present in the Taylor ex- 
pansion for any of the chosen h x allows us to introduce 
with sufficient accuracy the effective formula (JsJ) , which is 
desirable because of a smaller number of free parameters. 
The presence of two (nonconstant) terms in the Taylor 
expansion at small 1/s results in a noninteger exponent 
k in Eq. (JsJ) , where only one power is considered. The ex- 
ponent k is intermediate between the two corresponding 
integer exponents appearing in the Taylor expansion. 
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FIG. 2. The dependence of the magnetization (extrapo- 
lated to -/Vrg — > 00) on the inverse number of states kept, 1/s, 
at 7 = 1, in the absence of randomness, ho = 0. The solid 
lines are fits of the form M = M(0) + B x (l/s) fc , the values 
of the exponent k are shown. 

Finally, (v) we extrapolate (00, h x , ho) to zero ex- 

ternal field h x , using the expression M (po,h x ,ho) = 

M °°' ) (00, 0, ho) + Ch xl where C is a fitting parameter, 

Mo(ho) = M (00, 0, ho) is the spontaneous magnetiza- 
tion as a function of ho (at fixed 7) and the exponent a 
is determined by the fit to numerical data and is smaller 
than 1 in the Griffiths' region, where the susceptibility 
diverges as h x — > 0. The above formula is a slight mod- 
ification of Eq. (||), based on the observation that the 
exponent k is negligible. The modified formula repro- 
duces the behavior of numerical data for any value of the 
factor of anisotropy 7. As Fig. [3] shows, the magnetiza- 
tion is decreasing with increasing ho, and we clearly see 
that Mq goes to zero at 7 = 1 even though the uncer- 
tainty is rather large which is due to that we are exactly 
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at criticality, e.g. M (h = 0) = -0.18 ± 0.18, when 
properly extrapolated. 
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FIG. 3. The dependence of the magnetization (extrapo- 
lated to TVrg — > oo, s — > oo) on the uniform magnetic field 
h x at 7 = 1. The lines are fits of the form M — Mq + Ch". 

We point out that the above successive extrapolations 
are not needed in the Ising limit (7 = 0) Ref. [l] and add 
but little accuracy for small 7. This is due to the fact that 
the truncation error decreases very rapidly as a function 
of the number s of states kept when 7 is small. Thus, 
provided the system size is reasonably large and s not too 
small, the only extrapolating procedure needed to get the 
correct behavior for the spontaneous magnetization, close 
to the Ising limit, is h x — > O.El 

IV. RESULTS AND DISCUSSION 

As a starting test for the procedure discussed in Sec. 
Ill, we have calculated the spontaneous magnetization 
Mq for many different values of 7 at zero random held. 
As we anticipated, quantum fluctuations tend to reduce 
the magnetization. We fitted our data as 

M (/ l0 = 0,7) = i(l-7 2 ) ,? , (9) 

with d = 0.255t°;oJ5, very close to 1/4 (the circles 
in Fig. which is to be compared to the SW re- 
sult Eq. ([§ (the dashed line in Fig. |). In SW lan- 
guage, the factor (1 — ^ 2 Y represents the reduction of 
the magnetization, with respect to the saturation value 
Mo(ho = 0,7 = 0) = 1/2, due to quantum fluctuations. 
This factor is independent of the sign of 7, and drives the 
spontaneous magnetization to zero as 7 — > 1. Thus, even 
in the absence of randomness, a critical behavior is in- 
duced by quantum fluctuations. This behavior interferes 
with the critical behavior controlled by the random held. 
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FIG. 4. The triangles represent the critical points h c 
obtained by DMRG calculations for different values of the 
factor of anisotropy 7. The solid line is the best fit of the 
form h c = h(0)(l - -y 2 ) e with g = 0.52 ± 0.01. The system 
is in a ferromagnetic state below this line. The circles rep- 
resent the dependence of the spontaneous magnetization on 
7 at h = 0. Data are htted by M = 0.5(1 - 7 2 ) 1 ' with 
i9 = O^SlQ'oJg . The error bars that are smaller than the 
symbol size are not shown. The result of SW theory, Eq. (Q) , 
is also drawn (dashed line). 

To analyze this interplay, we applied the same pro- 
cedure to obtain the spontaneous magnetization in the 
presence of randomness (Fig. ||). Only s = 8 states 
proved to be enough to obtain accurate spontaneous- 
magnetization data for 7 = 0.4 and no additional extrap- 
olation to s — > 00 was necessary. The above mentioned 
procedure was fundamental, instead, to produce sensible 
data for 7 = 0.8. About 1000 runs were needed for a 
given strength of the random held ho, the uniform field 
h x and the number of states kept s to get an accurate 
statistical average of the magnetization. The available 
computational resources limited us to s < 12, so it was 
very important to have the fit formula (^|) with a small 
number of free parameters. The accuracy of the final Mq 
data is good despite the small number of states s used 
in the calculations, which induces errors in the exponent 
k and consequently in M(l/s — 0), (see Eq. and the 
text below it). 

Once the spontaneous magnetization is obtained, the 
next step is to determine the phase-transition line h c ("f), 
i.e. the line in the ho vs 7 plane where the spontaneous 
magnetization M vanishes. Since the data close to criti- 
cality are usually affected by large errors, the best way to 
obtain h c is to find Mq in the ordered phase and then fit 
the data to formula <JsJ) , leaving j3 ~ 2—tfi as an adjustable 
parameter: Mo = Dx (h c — ho) 13 , as it is shown in Fig. (||). 
The results are h c = 1.234 ± 0.007 and f3 = 0.38 ± 0.03 
for 7 = 0.4 and h c = 0.79 ± 0.01 and (3 = 0.35 ± 0.04 for 
7 = 0.8. We emphasize here the fact that the magnetiza- 
tion fit gives an exponent (3 compatible with 2 — </> « 0.38 
for any 7, which is a strong numerical evidence that the 
system falls in the same universality class as the quan- 
tum Ising spin chain for all 7 < 1. In the language of 
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SW theory the presence of a gap in the excitation spec- 
trum leads to an effective purely Ising-like model, with a 
renormalized coupling constant J 7 , near criticality. 
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FIG. 5. The dependence of the spontaneous magnetiza- 
tion on the strength of the random magnetic field ho- The 
data points were fitted by the formula Mo = D x (h c — h{ 
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The phase diagram is plotted in Fig. (||). The critical 
line (the one passing through the triangle symbols which 
correspond to our numerical results) is well approximated 
by the equation 



M7) = M0)(i- 7 2 ) £ 



(10) 



where h c (0) = e/2 and g = 0.52 ±0.01, very close to 1/2. 
Thus the effect of quantum fluctuations in the Griffiths' 
region around the critical point is such as to increase the 
effective strength of the random field. In other words, as 
the magnetic order is weakened by quantum fluctuations, 
a weaker random field ho is needed to drive the system 
to the paramagnetic phase. We point out that the crit- 
ical field strength vanishes for 7 — > 1, where a different 
critical behavior sets in [cf. Eq. @], controlled by the 
corresponding critical point. 

In the language of SW theory quantum fluctuations 
lead to a renormalized effective Hamiltonian. The re- 
sult (^|) indicates that the field Sf is renormalized as 
Sf = (1 - ^fSf at ho = 0. This in turn leads to an 
effective Ising-like model with renormalized coupling con- 
stant J 7 = J(l — 7 2 ) 215 . Thus the result ( |i~0| ) indicates 
that the increase in the strength of random field due to 
quantum fluctuations is entirely due to a reduction of 
the coupling constant, while the coupling to the random 
field is not renormalized, so that the scaling law g = 2i9 
holds. In other words, the interference of the two critical 
behaviors leads to the relation (2Mo) 2 = 2h c /e, i.e. the 
ratio Mg/hc = l/2e is independent of 7. 

To check the consistency of the phase diagram we have 
calculated the energy gap G between the ground state 
and the first excited state, which provides an indepen- 
dent determination of the critical point. Indeed, ran- 
domness tends to fill in the gap, which should become 
zero at criticality. This is a property, which is easier to 



calculate, compared to the spontaneous magnetization, 
because it can be calculated at zero uniform magnetic 
field and no extrapolation for h x —> is needed. How- 
ever, unlike the ground-state properties, the properties 
of the excited states require the full extrapolating pro- 
cedure -/Vrg — ► 00, s — ► 00, to be properly determined, 
even in the Ising limit. Indeed the energy gap depends 
strongly on the number s of states kept in the DMRG 
procedure for all 7. Some examples are shown in Fig. ||, 
in the absence of randomness. It is seen that, once the 
correct extrapolation procedure is adopted, the DMRG 
reproduces the exact result (^) within good accuracy. 
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FIG. 6. The dependence of the energy gap on the in- 
verse number of states kept in the DMRG truncating proce- 
dure. The system size is up to 404 sites, h x = 0, ho = 0. 
The straight lines are fits of the form G = a + b/s, where 
a = 0.605 ± 0.003 for 7 = 0.4 and a = 0.197 ± 0.004 for 
7 = 0.8. 

We have calculated the energy gap G for different 7 
and ho (Fig. 7j). As ho increases, G is reduced from its 
initial value ( 3) and is driven to zero as the critical point 
is approached. The results for the gap are in agreement 
with the results obtained from the spontaneous magne- 
tization data (cf. the triangles in the Fig. ||). The curves 
are exponentially flat near criticality so that the crit- 
ical points are much more accurately determined from 
the magnetization curves. However, once h c is deter- 
mined from magnetization data, the data corresponding 
to — 1/ ln(G/Go) collapse near criticality on one curve as 
a function of 1— ho /h c , as shown in the inset of the Fig. [7]. 



G 




FIG. 7. The dependence of the energy gap on the strength 
of the random field ho- The system size is above 200 sites, 
h x = 0. The average is performed over 500 realizations of the 
random field. The inset shows rescaled data, and the solid 
curve is -l/ln(G?/G ) = K(l - h /h c ) x with K = 2.6 ± 0.4 
and A = 1.6 ± 0.2. Go denotes the gap in the absence of 
randomness, Eq. (ji|). 

So far, we have discussed the procedure to obtain the 
critical behavior of the anisotropic XY model in a ran- 
dom field, and have studied those quantities which vanish 
at criticality. Now we analyze the correlation length and 
we look for the nearly critical behavior at a (large) fixed 
system size and a given number of states in the truncation 
procedure kept. In particular, since there is a great deal 
of analysis on finite-size effects in approximate numerical 
calculations for models near criticality, we concentrate 
on the distinguishing feature of the DMRG method, i.e. 
the effect of the truncation of the Hilbcrt space. Thus 
we calculated both the connected and nonconnected two- 
point spin-spin correlation functions along the easy axis 
of magnetization 



Cj(r) = 4[(S ^; 
C,-(r) = 4<SJS?>, 



(11) 

(12) 



where ( • ) stands for the quantum expectation value in 
the ground state, for any given realization j of the ran- 
dom field, and we are using here an abbreviated notation 
with respect to Eq. ([?]). A factor of 4 in the definition 
of Cj (r) was introduced to allow for a direct comparison 
with Ref. in the Ising limit 7 = 0. We point out that, 
for a given realization of the random field, (Sq ) and (Sy ) 
are different, and neither of them is equal to the magne- 
tization given by Eq. (Q), due to the lack of translational 
invariance. 

We used a method which is slightly different from (but 
essentially equivalent to) the method described, e.g., in 
Ref. [l]. We grew the system up to a size of ~ 400 sites, 
while storing the truncated spin operators for the differ- 
ent sites of the chain. The typical Cj(r) ^should decay 
exponentially with a correlation length £.13 The plot of 
the average over disorder A" 1 *£j=i In Cj(r) as a func- 
tion of the distance r allows then to extract the harmonic 




FIG. 8. The average (over 500 configurations) of the log- 
arithm of the connected correlation function for 7 = 0.8. The 
typical correlation length £ is calculated as the negative in- 
verse slope of these curves. 

average of the correlation length from the slope of the 
linear behavior at large r (see, e.g., Fig. . By plotting 
the £ vs ho for each 7 (Fig. m, we obtained the peaks that 
indicate the would-be critical point. As it is evident, the 
peaks are rather broad, and the correlation lengths rather 
short, due to the truncation of the Hilbert space, but the 
position of the maximum can be extracted with a good 
accuracy, and we found /i pGa k = 1-40 ± 0.05 for 7 = 0.0, 
Vak = 1-25 ± 0.05 for 7 = 0.4, and /i pcak = 1.00 ± 0.05 
for 7 — 0.8. We checked that increasing the number of 
states s kept in the truncating procedure leads to nar- 
rower peaks which shift towards the critical values h c (^f) 
that were obtained from the spontaneous magnetization 
data in the infinite-size and s — > 00 limit. 
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FIG. 9. The dependence of the typical correlation length £ 
on the strength of the random magnetic field ho- There were 
s — 8 states kept in the truncation procedure for the white 
data points. 

To see how close to criticality the system is at the peaks 
of the correlation lengths we looked for a typically crit- 
ical property of the correlation functions. In Fig. ^ we 
plotted the average over the realizations of the random 
field of the logarithm of the correlation function (|lT 
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In C(r) = iV c 1 J2j=i m< 7j'( r ) as a function of y/r, since 
such a behavior is expected at criticalityBB Our result 
shows that the curves are very sensitive to the value of 
ho, and bend upwards when ho is smaller than some par- 
ticular value h or downwards when ho > h. The best 
fit to linear behavior is found at h = 1.44 for 7 = 0.0, 
h = 1.30 for 7 = 0.4, and h= 1.03 for 7 = 0.8. These 
values and the corresponding ft. pca kj obtained from the 
£ vs ho curves, coincide within the error bars. Thus the 
truncated Hilbert space embodies the critical behavior of 
the system in a self-consistent though approximate way. 
Increasing the size of the Hilbert space improves the ac- 
curacy in the description of the critical properties until 
the only limitation to a fully developed criticality is the 
finite size of the system. It is perhaps important to re- 
mark that the extrapolating procedure proposed in Sec. 
Ill to extract the spontaneous magnetization may also be 
applied to the correlation functions. The only difference 
is that much more computational time and memory is re- 
quired to calculate and store the data which refer, for a 
given distance r, to different system sizes, different num- 
ber of states s, and different random field realizations. 




FIG. 10. The average of the logarithm of the noncon- 
nected correlation function. The curves correspond to the 
following strengths of the random magnetic field (from top to 
bottom): h = 1.25; 1.44; 1.65 for 7 = 0.0 (solid lines), 
ho = 1.1; 1.3; 1.5 for 7 = 0.4 (dot-dashed lines), and 
h = 0.8; 1.03; 1.2 for 7 = 0.8 (dashed lines). 

V. CONCLUSIONS 

In summary, we analyzed some properties of the spin- 
1/2 quantum anisotropic XY chain in a transverse ran- 
dom magnetic field by means of the density-matrix renor- 
malization group. The dependence of the magnetization 
on the uniform magnetic field h x , the strength of the 
random magnetic field ho and the factor of anisotropy 7 
was obtained. The order-disorder phase-transition line 
was determined [Eq. (|l0|)] and the phase diagram was 
drawn. The energy gap between the ground state and 
the first excited state was investigated as well. The de- 



pendence of the gap on both the strength of the random 
field ho and the factor of anisotropy was obtained. The 
gap vanishes at the phase transition, determined inde- 
pendently from magnetization data, and reproduces the 
correct limiting value (||) as the strength of the random 
field is reduced. Finally we calculated the connected and 
nonconnected two-point spin-spin correlation functions 
along the easy axis of magnetization. We studied in par- 
ticular the asymptotic behavior for large distances at critg 
icality, and found that the -y/r-behavior found at 7 = OU'u 
persists for 7 > 0. 

The critical properties are remarkably the same for all 
7 < 1, providing clear numerical evidence for univer- 
sality, i.e. models with different values of the factor of 
anisotropy < 7 < 1 all belong to the same universality 
class as the spin-1/2 quantum Ising chain in a transverse 
random magnetic field. 

The main advantages of the method used in this pa- 
per to investigate the properties of a random quantum 
system, with respect to other numerical methods have 
been discussed in our previous paper.EJ Here we wish to 
comment in deeper detail on the specific technical prob- 
lem which was dealt with in this paper. Indeed, in the 
present case, the density-matrix renormalization-group 
approach requires a particular care, due to the interplay 
of randomness, finite-size effects, and because the conse- 
quences of truncation of the Hilbert space are enhanced 
by the presence of quantum fluctuations (with respect to 
an Ising-like ground state) as 7 is increased. This effect is 
most easily seen in the gradual spread of the eigenvalues 
of the density matrix, i.e. in the increasing importance 
of including more states in the truncating procedure, in 
order to obtain an accurate description of the system. 
However a reliable protocol was discussed in this paper, 
which takes care of all those aspects on equal footing, al- 
lowing us to obtain physical results which are very robust 
with respect to "local" variations of the extrapolating 
procedures. 

Criticism was often raised against the use of the 
DMRG in the investigations of random systems (see, e.g. 

and some modifications of the DMRG procedure 
were proposed to deal with randomness. c3 The authors 
essentially refer to the failure in accommodating sudden 
changes of the ground state within a truncated basis for 
the Hilbert space. These objections are appropriate, in 
principle, and suggest a careful analysis of the stability of 
the DMRG results as the size of the basis is enlarged. In 
this paper we showed that this analysis is possible, and 
the objections may be overcome. 
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